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The vibrational branching ratios and asym- 
metry parameters for C0 2 have been 
determined in the wavelength region of 
650 A to near the ionization onset at 
about 840 A. The study was performed 
using synchrotron radiation from the 
Daresbury storage ring that was dispersed 
with a 5 m grating monochomator that 
afforded resolution of 0.1 A to 0.2 A. This 
resolution allowed the study of the 
branching ratios and asymmetry parameters 
with enough detail to see the changes in 
the parameters within the pronounced 
autoionization structure in C0 2 in this 
wavelength region. While the electron 
spectrometer resolution was not sufficient 
to resolve the spin orbit and Renner- Teller 
splitting in the photoelectron spectra, we 
are able to fit the data with a model that 



identifies the major structure in terms of 
the symmetric stretch and elements of the 
asymmetric stretch and bending modes. 
A calculation of the expected relative 
vibrational excitations based upon the 
Franck-Condon principle clearly showed 
non-Franck-Condon behavior in some of 
the vibrational-electronic transitions. 
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Introduction 



The study of C0 2 photoionization has historically 
attracted much attention due to the importance of the 
ionization of C0 2 in the photophysics of planetary 
atmospheres, including the Earth's atmosphere [1,2]. 



Additionally, C0 2 is an integral part of the carbon cycle 
for plant life and, as a consequence, the photochemistry 
and photophysics of this molecule has been the subject 
of considerable interest by a number of scientists over 
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the last several decades. The ready availability of high 
purity samples of C0 2 and its ease of use in various 
experimental arrangements have produced many experi- 
mental efforts studying its various properties [3-7]. 

The ionization onset region in the photoabsorption of 
C0 2 contains a number of strong autoionization features 
that converge on the A 2 II U and B 2 2* states seen in a 
photoionization mass spectrometric study by McCulloh 
[8] who tabulated them and gave the various series 
identifying labels which we will use in this work. 
Hubin-Franskin et al. [9] and Holland and Hayes [10] 
have used photoionization techniques coupled with 
synchrotron radiation sources to study the branching 
ratios or vibrational transitions in the wavelength region 
reported in this study [9,10]. Baer and Guyon [11] 
used synchrotron radiation to study autoionization 
and isotope effects in the photoionization in the 620 A 
to 900 Awavelength region. These researchers used a 
threshold photoelectron spectrometer system and identi- 
fied a number of vibrational transitions that are ordinar- 
ily weak or forbidden by selection rules. 

Recently Buenker et al. [12] have performed calcula- 
tions and reviewed the status of the studies on C0 2 
electron energy loss spectroscopy. This technique can 
give results for cross sections and oscillator strengths 
that are similar to those obtained by photoionization and 
can be a direct test of calculational techniques. Takeshita 
et al. [13] recently calculated the ionization energies and 
the Franck-Condon factors for C0 2 photoionization. 
Johnson and Rostas [14] reviewed the spectroscopic 
literature for the vibronic structure of the ground state 
and first few excited electronic states of C0 2 . The 
spectroscopic results given by these authors were most 
useful in developing a method of analyzing the data 
reported in this work. 

In a previous publication [15], we gave results for the 
branching ratios and asymmetry parameters for C0 2 in 
the wavelength region of 685 Ato 795 A for some of the 
autoionizing lines in the Tanaka-Ogawa (TO) series, the 
Henning sharp (s) and diffuse (d) series. However, we 
did not publish the extensive additional data we had 
obtained as a result of analyzing the photoelectron 
spectra taken in this wavelength region. The results 
reported here are a complete analysis of all our data, 
providing extensive information on the photoionization 
behavior of C0 2 in this wavelength region. Another 
publication from our collaboration used the angular 
distribution information to determine the electronic 
structure, symmetry, and decay dynamics for members 
of the TO series [16]. It was noted in these publications, 
as well as in preliminary work done in the early 1980s, 
that the strength of the vibrational progressions were not 
described by a simple Franck-Condon model of photo- 
ionization [17]. We also explored the effects of vibronic 
coupling on the intensity distribution in the 4 cr g 



photoionization channel of C0 2 [18]. Hubin-Franskin et 
al. studied C0 2 autoionization and pointed out that the 
use of the method worked out by Smith and modified by 
Eland to apply the Fano-Mies configuration interaction 
theory to molecules could be used to explain portions of 
the non-Franck-Condon behavior seen in the photo- 
ionization of C0 2 [19-23]. These papers pointed out the 
difficulty of completely accounting for the vibrational- 
electronic transition intensities due to configuration 
mixing of electronic states, vibronic and rotational inter- 
actions, and the competition between the many paths 
available for an excited state neutral molecule to decay. 
The Smith formalism requires a knowledge of the 
Franck-Condon factors for all the transitions involved, 
ground state to autoionizing level and autoionizing level 
to all the possible final states. While this calculation is 
possible with diatomic molecules, it is much more 
difficult with larger molecules with additional degrees 
of freedom. In the case of C0 2 , the situation is even 
more complicated due to Renner-Teller splitting of the 
electronic levels when the bending mode is excited. 
As a result, there has been no systematic theoretical 
exploration of the effects of autoionization on the 
vibrational branching ratios and asymmetry parameters 
forC0 2 . 

In our previous publication concerning the branching 
ratios and asymmetry parameters in this wavelength 
region we presented only a small subset of the data we 
have on C0 2 in this wavelength region. We present here 
a detailed explanation of the data accumulation 
techniques and a complete discussion of the analysis of 
the data that was not possible in the shorter summary of 
this work [15]. To be consistent with the earlier NIST 
publications on this topic, we have used the unit of the 
Angstrom to describe the wavelength scale used for 
these experiments even though this unit is no longer 
sanctioned as an acceptable metric system unit. 

2. Experimental Procedure 

This experiment was carried out on the 5 m normal 
incidence monochromator fitted to a beamline at the 
Daresbury SRS, providing a photon flux of ~10 10 pho- 
tons/s within a bandpass of 0. 1 A in the spectral region 
from 650 A to 840 A [24]. The light was brought into 
the experimental chamber by a 2 mm internal diameter 
glass capillary light guide whose capillary aperture was 
placed close to the exit slit of the monochromator. In 
addition to offering a low-loss transport for the vacuum 
ultraviolet radiation, the capillary also served to main- 
tain a pressure differential between the experimental 
chamber and the ultra high vacuum of the optical 
monochromator. A schematic diagram of the experi- 
mental apparatus is shown in Fig. 1. The light guide 
extends from the exit slit, EX in Fig. 1, to the interaction 
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Fig. 1. Cross section view of the electron spectrometer chamber. The chamber is 76 cm in diameter and is 91.4 cm in length. The major aspects 
of the instrument are labelled in the figure. 



region above the gas entrance tube, GS in Fig. 1, and can 
be as long as 30 cm dependent upon exact placement 
of the experimental chamber with respect to the mono- 
chromator. 

The electron spectrometer system comprises two 
100 mm mean radius hemispherical analysers, one 
rotatable about the incoming light beam as an axis and 
the other fixed, contained in a chamber shielded from 
magnetic fields by three layers of |x-metal. The radiation 
from the monochromator is polarized, its polarization 
depending upon the optics of the monochromator and 
the subtended angles of acceptance of the synchrotron 
radiation. In our configuration, the light is polarized 
with between 75 % and 80 % of the light intensity 
having its electric field vector (E-vector) perpendicular 
to the plane of the paper in Fig. 1. In this configuration, 
the fixed analyser accepts electrons ejected parallel to 
the E-vector of the incident radiation. The fixed analyser 
is ES-2 in Fig. 1. The other electron spectrometer is 
rotatable about an axis defined by the direction of the 
light and hence collects electron in the plane of the 
E-vector. This movement allows the angular distribution 
of the photoelectrons to be explored completely and is 



sufficient to determine the angular asymmetry para- 
meter for the scattering process. The entrance lenses for 
the spectrometers are three element zoom lenses based 
upon the design of Harting and Read [25]. The entrance 
cone to the lens system has a small aperture, usually 
about 1 mm in our experiments, which acts as the 
limiting aperture for determining both the energy and 
angular resolution of the system. The zoom lens focuses 
the electrons from a small interaction volume deter- 
mined by the size of the light beam exiting the capillary 
and the size of the gas jet exiting the gas entrance tube 
onto the entrance plane of the hemispheres. The pass 
energy of the electron analyzer and focus voltages are 
set by external controls. The pass energy remains fixed 
for a particular set of experiments, and the other 
voltages are appropriately varied to scan the electron 
energy spectrum as required by using an automated data 
control system. The electrons are dispersed upon 
passing through the analyzer hemispheres and focussed 
on the hemisphere exit plane. Since the apparatus was 
first described in the literature, the electron spectro- 
meters have been fitted with position sensitive detectors 
which are placed near the exit focal plane of the 
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hemispheres [26]. This allows for the simultaneous 
detection of a range of energies in the photoelectron 
spectra and thereby improves the data quality for a given 
period of data accumulation, compared to using a 
conventional electron multiplier behind an exit slit. 

The polarization of the incoming light was measured 
using a three mirror polarizer with tungsten mesh and 
plate photodiodes which could be rotated with the 
rotatable analyzer through 90° in order to determine the 
light polarization. The polarization detection device was 
constructed based upon considerations given by Horton 
et al. [27]. The polarization was checked frequently 
since small movements in the storage ring beam 
position, and the mirrors focussing the light onto the 
entrance slit of the monochromator, can have a marked 
effect on the polarization. It was found that, provided 
the pre-mirror adjustments were kept optimized for 
maximum photon flux at the exit slit of the mono- 
chromator, the polarization would remain stable during 
the accumulation of a particular data set. The tungsten 
wire mesh at the entrance of the polarizer served as a 
photocathode for monitoring the intensity of the incom- 
ing light beam. A tungsten plate serving as a second 
photocathode collected the beam remaining after three 
reflections of the beam. The ratios of these photo- 
currents as a function of the angular position of the 
movable electron spectrometer provided the data 
necessary to determine the polarization of the light. 

Calibration of the energy response of the analyzers 
was performed using the known values of the cross 
section and asymmetry parameters for argon or helium 
gas and following standard procedures outlined in the 
literature [28-30]. For all the spectra reported here the 
electron spectrometer resolution was determined from 
the rare gas calibration to be 4 1 meV for the fixed analy- 
ser and 46 meV for the rotatable analyzer. The 5 m 
monochromator resolution was 0.1 A(~2 meV) for the 
measurements taken at wavelengths shorter than 750 A. 
At wavelengths longer than this, where the structure in 
the absorption spectrum is less dense, the resolution 
requirements could be relaxed and a wavelength resolu- 
tion of 0.2 A was used. The data were accumulated 
by simultaneously taking photoelectron spectra at two 
angles with respect to the polarization direction by 
utilizing both electron spectrometer systems. Since the 
two analyzers could be positioned at different angles, 
a particular data point did not require rotation of the 
movable electron spectrometer system. The wavelength 
on the monochromator was then incremented by 0.1 A 
and another set of photoelectron spectra taken. The light 
polarization was checked periodically by a 90° rotation 
of the movable electron spectrometer. During data 
accumulation, the time spent at a particular electron 
kinetic energy was determined by integration of the light 



flux signal to some predetermined amount so that all 
the points in a particular data set would be correctly 
normalized to the same total light flux. 

The differential cross section for photoabsorption in 
the dipole approximation for a randomly oriented gas 
may be expressed as 



d<r„ cr v B v 

d72 = 4^ [1+ T (3/,cos20+1) 



(1) 



where 6 is the angle between the major polarization axis 
and the ejected electron, P is the degree of polarization 
of the incoming light, (I is the solid angle of collection 
of the photoelectrons, and a v is the partial cross section 
for the vibrational-electronic channel corresponding to 
the photoelectron being detected. The total cross section 
for a particular electronic channel is the sum of the 
partial cross sections of the individually resolved 
vibrational channels. The vibrational branching ratio is 
defined as the partial cross section for that 
channel divided by the total cross section for the elec- 
tronic state. The number of photoelectrons per steradian 
per unit light flux, dNJ&fl is proportional to the differ- 
ential cross section and hence we can recast the above 
equation into one that refers to the measurable para- 
meters of the experiment: 



&N V N v B v 

-T7T = -r- [1+^r (3P cos 20+1) 

d/2 4ir 4 



(2) 



N v is the total number of electrons detected in a partic- 
ular vibrational channel summed over the all solid 
angles. Measurements were made simultaneously at 
= 0° and 6 = 90° and hence B„ and N„ could be directly 
deduced from the two spectra [31]. The branching ratio 
for a particular transition is the ratio of N v with respect 
to the sum of all the N v 's for a particular electronic 
excitation. The data reported here were taken only for 
transitions that left the COJ molecule in the X 2 II g 
ground electronic state. Altogether about 1500 data sets 
were taken in the wavelength region 650 A to 890 A. 

Figure 2 shows the photoionization efficiency (rela- 
tive photoionization cross section) for CO! in the wave- 
length region of interest in the present study. The data 
were taken using a rotationally cooled sample 
of C0 2 and a laboratory light source coupled to a 
quadrupole mass filter [16,32]. A 1 meter near-normal 
incidence monochromator provided dispersed radiation 
with a wavelength resolution of 0.12 A. This rotation- 
ally-cooled spectrum shows more detailed and sharper 
structure, particularly near ionization onset, than does 
the spectrum obtained by Berkowitz [33] with a slightly 
better wavelength resolution of 0.07 A. The spectrum 
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Fig. 2. Photoionization efficiency spectrum for CO2 in the wavelength region from ionization onset to 600 A,. The data were taken using a 
photoionization mass spectrometer system with a wavelength resolution of 0.12 A. 



shows members of the TO series that have as limits the 
A state of COJ vibrational levels. The notation An,, (TO) 
n = 4,5, . . etc., means the level is a member of the TO 
series having the n v vibrational level of the A state of 
COJ as its limit. The symbol B« (s,d) n = 3,4," etc. is a 
Rydberg level of quantum number n having as its series 
limit the B state of COJ in the vibrational ground state. 
The s or d refers to the Hennings sharp or diffuse series. 
The notation A3„(L) v = 0, . . refers to a level of the 
Lindholm series with principal quantum number 3 
that has the A state as its limit and has a vibrational 
excitation of v in the symmetric stretch mode. 

3. Analysis of the Data 

The analysis of the large amount of data taken on C0 2 
required careful organization of the data and the devel- 
opment of a reliable method of extracting the highly 
detailed information on vibrationally resolved asym- 
metry parameters and branching ratios that the data 
contained. The raw data must first be corrected for the 
differing detection efficiencies of the two electron 
energy analyzers. This correction is a function of the 



kinetic energy of the electron and is accounted for by the 
determination of an energy dependent correction factor 
from angular distribution data measured by the two 
electron spectrometers for rare gases as discussed in the 
previous section. The correction function was frequently 
checked during the course of the experiment by measur- 
ing the photoelectron spectra of a rare gas as a function 
of wavelength over a wavelength range sufficient to 
generate photoelectrons with kinetic energies equivalent 
to those produced in the photoionization of C0 2 . The 
system exhibited excellent stability in this regard as, 
once the lens voltages and other operational parameters 
of the electron spectrometers were set, the relative 
efficiencies of the analyzers did not vary in any 
significant manner. 

The 2 ll g ground state of COJ is spin-orbit split into 
two levels that are 19 meV apart and the bending mode 
populated in the resonance regions is split a similar 
amount by Renner-Teller coupling. This structure was 
not resolved and hence the effect of these interactions 
was to broaden the spectral peaks observed in the COJ 
photoelectron spectra. This could be accounted for by 
examining the shape of the photoelectron peaks seen in 
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the rare gas spectra and appropriately adding two such 
peaks that are offset by the 19 meV splitting causing the 
broadening. Using this procedure and some testing on 
fitting some actual COJ spectra resulted in obtaining 
47 meV and 56 meV as the best effective resolutions of 
the two analyzers to be used in the fitting routine. The 
fitting routine allowed for the resolution to be a free 
parameter or a constrained parameter but the best over- 
all fittings were obtained when the resolutions were 
fixed at these values and not allowed to vary during the 
fitting process. The fitting routine used as a model a 
series of peaks whose widths were fixed and whose 
relative positions were fixed by the spectroscopic 
positions of the vibrational energy levels. The fitting 
routine was a nonlinear least squares analysis [34] 
that was customized at NIST for the type of spectra 
observed. 

The spectroscopic values of the vibrational para- 
meters were taken from the literature [35,36] and fitted 
to our data using &>i = 156.7 meV for the symmetric 
stretch, co 2 = 0.060 meV for the bending mode, and 
co 3 = 182.3 meV for the asymmetric stretch. These 
values gave good fits to the data over a broad range of 
electron kinetic energies and represent an averaging over 
the levels split by the Renner-Teller interaction. The 
transitions to the asymmetric stretch and bending mode 
with odd number of quanta are forbidden by 
selection rules from the ground state of the neutral 
molecule but have been observed in other photoelectron 
work reported in the literature [37]. Though these 
transitions are not allowed by the Franck-Condon 
principle, their existence is well established experimen- 
tally and can be attributed to anisotropic interactions of 
the escaping electron with the molecular core. The 
presence of these transitions and the Renner-Teller split- 
ting of the levels gives rise to an overlapping of the 
possible transitions, and this makes unique identification 
of the photoelectron peaks we observe impossible ex- 
cept in some prominent circumstances. Chambaud, 
Gabriel, Rosmus, and Rostas have identified many 
of the low lying levels of COJ, providing a basis for 
estimating the positions and compositions of the levels 
we observe [38]. Using the notation (vi,v 2 ,v 3 ) where v\ 
is the quantum number of the symmetric stretch, v 2 
the quantum number of the bending mode, and v 3 that 
of the asymmetric stretch, we find the spectral features 
can be adequately assigned to the following vibrational 
levels: 

(vi,0,0) vi, = 0, 1... 
(vi,l,0) vi, = 0, 1... 
(vi,2,0) vi, = 0, 1, ...4 
(vi,0,l) vi, = 0, 1, ...5 
(v,,l,2) Vi, > 3 



The upper value of Vi is determined by the span of the 
energy range covered by the photoelectron spectra. In 
some cases it was not necessary to include the (vi,0,l) 
or (v b l,2) levels, as the intensity in the high binding 
region was small and the structure was adequately repre- 
sented by the simpler set of levels. Table 1 gives the 
vibrational levels used in the fitting model in the energy 
range of one eV of vibrational excitation. The levels are 
identified as discussed above and the energy position is 
given with respect to the (0,0,0) vibrational ground state 
of COJ. 

Table 1. The vibrational levels used in the least squares fit process 



Vibrational level 


Energy (eV) 


Vibrational level 


Energy (eV) 


000 


0.0000 


420 


0.7474 


010 


0.0603 


500 


0.7835 


020 


0.1206 


401 


0.8091 


100 


0.1567 


510 


0.8438 


001 


0.1823 


411 


0.8691 


110 


0.2170 


312 


0.8944 


120 


0.2773 


600 


0.9402 


200 


0.3134 


501 


0.9658 


101 


0.3390 


610 


1.0005 


210 


0.3737 


511 


1.0261 


220 


0.4340 


412 


1.0511 


300 


0.4701 


700 


1.0969 


201 


0.4957 


710 


1.1572 


310 


0.5304 


512 


1.2078 


211 


0.5560 


800 


1.2536 


320 


0.5907 


810 


1.3139 


400 


0.6368 


612 


1.3645 


301 


0.6524 


900 


1.4103 


410 


0.6871 


910 


1.4706 


311 


0.7127 


712 


1.5212 



It is possible other types of excitations may contribute 
to the spectra observed. The appearance of finite values 
for intensity at an energy determined by the above en- 
ergy levels does not constitute an unambiguous identifi- 
cation of a particular vibrational level for COJ. As an 
example, our data would not separate out the (2,0,0) 
components from those of components of the split 
(1,0,1) and the (1,2,0) levels. The (1,2,0) level is split 
into 4 levels by the Renner-Teller and spin-orbit interac- 
tions, the upper two of which are approximately coinci- 
dent with the (2,0,0) levels. The electron energy resolu- 
tion of the present experiment was not sufficient to 
warrant attempts at resolving the levels due to the 
splittings. Similarly, the intensity observed for 
"forbidden" transitions to levels of type (n,l,0) should 
not be regarded as definitive. In an earlier publication- 
where we studied C0 2 photoionization with an electron 
spectrometer resolution of about 19 meV, members of 
the (vi,0,2) progression were clearly visible [16] which 
indicates that the asymmetric stretch excitations can 
play a role in the photoabsorption of C0 2 . In this lower 
resolution study, these excitations could not be resolved 
from those excitations shown in Table 1. Our goal in 
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the fitting was to identify peak positions which would 
best account for the electron intensity observed and 
which could be identified with predicted vibrational 
structure. 

To accomplish the fitting, relative peak positions were 
fixed by the spectroscopic model chosen, the peak 
widths were fixed as determined in the previous discus- 
sion and the heights of the photoelectron peaks were 
allowed to vary in order to minimize the fitting errors 
with the least squares technique. After many trials it 
soon became evident which vibrational modes were the 
most important to achieve a reasonable fit. This set of 
vibrational modes was adopted as a basis for the analy- 
sis. At energies 2 eV above the threshold, it was not 
possible to resolve the individual vibrational peaks and 



they were often seen as a modulated background whose 
intensity was accounted for by high members of the 
above progressions. In the wavelength region of 
652.08 A to 680.98 A it was sufficient to use only the 
first three vibrational progressions shown above; in the 
region 689.7 A to 713.9 A it was necessary to use the 
first four levels; and in the remainder the complete set 
was used. The only exception to this was the use of 
levels of the sort (vi,l,l) in the 689.77 A to 694.64 A 
region where this addition seemed to improve the fit 
somewhat but did not suggest itself elsewhere and hence 
will not be pursued in this report. 

After fitting, each data set and the fit were inspected, 
and the model varied by adding or removing vibrational 
modes as required to achieve the best fit. Figure 3 
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Fig. 3. Two photoelectron spectra taken on the 0° analyzer at nearby wavelengths. The 
spectrum differ considerably but the same model fits both relatively well. The solid line 
is the fit and the dots are the actual data. The vertical columns represent the intensities of 
the various vibrational transitions as defined in the figure that contribute to the intensity 
at that point. 
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shows the quality of the fit obtained by the 0° analyzer 
at a photon energy of 17.180 eV. The lower portion of 
Fig. 3 shows the fitting results for the same analyzer at 
a photon energy of 17.218 eV. Both fits were accom- 
plished using the same model for vibrational state 
composition even though the structure is substantially 
different in the two closely spaced regions of the 
absorption spectrum. The vertical lines in the figures 
are the calculated amplitude of the contributions from 
the progression listed. The solid line, labelled FIT, 
represents the fit of the calculated spectrum to the raw 
data (DATA in the figure). The data points are for the 
most part obscured by the fitted line, an indication of the 
quality of the fit. 

It can be seen that the data are of sufficient quality to 
give reliable intensities for the vibrational components 
marked, allowing vibrational branching ratios down to 
the 5 % level to be measured. In the figures shown later, 
the uncertainties include the statistical element (type A 
uncertainty) derived from the fitting procedure, and 
further contributions due to uncertainty in the polari- 
zation of the incoming light and the calibration of the 
tungsten photodiodes (type B uncertainty) [39]. In 



general, the statistical uncertainty gave the largest effect. 
The error bars shown on the data represent the 
combination in quadrature of the uncertainties of 
type A and type B propagated with the appropriate 
error propagation for the equations given earlier in this 
paper. 

Having obtained the intensities for the vibrational 
levels for each of the two analyzers, we then used these 
data, together with the measured light polarization and 
the electron spectrometers' energy efficiency curves, 
to calculate the asymmetry parameters and branching 
ratios. The total electron count as a function of wave- 
length is shown in Fig. 4. This was obtained by sum- 
ming over the values of the calculated N„ for all the 
levels calculated for the X ground state. The intensities 
were corrected for the varying efficiency of the tungsten 
photodiode used to measure the light intensity and can 
be seen to agree closely with the results shown for the 
photoionization efficiency in Fig. 2 [40]. 

The three dimensional plot shown in Fig. 5 shows the 
electron spectra where the most intense members have 
been truncated in order to reveal the detail in the higher 
vibrational modes. Examination of this figure reveals 
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Fig. 5. Three dimensional plot of the results of the fitting routine for all the spectra measured. The molecular energy axis represents the vibrational 
excitation of the ionic ground state of COJ. The plot dramatically shows the increased vibrational excitation caused by autoionization phenomena. 



the amount of COJ vibrational excitation occurring as 
a result of autoionization which populate molecular 
energy levels well above 0.5 eV. This excess population 
of excited state CO! would have consequences on the 
reactivity of this molecule in atmospheres where COJis 
formed by uv irradiation. When the number of electrons 
as a function of molecular energy are summed as shown 
in Fig. 5, the total electron count as shown in Fig. 4 is 
obtained. 

4. Discussion 

In Figs. 6-20, the spectra for j3 and branching ratio are 
shown, for all the vibrational levels for which there was 
significant intensity. These have been selected from the 
complete data set because they show structure identi- 
fiable with the resonance positions; the complete data 
set is available from a NIST database as described 
below. 

To proceed with further analysis, we consider the 
Born-Oppenheimer approximation and the Franck- 
Condon (FC) principle. These two approximations have 
been used successfully for the analysis of vibrationally- 



resolved photoelectron spectra. Within these well- 
known approximations, the branching ratio in the 
photo-electron spectrum, for ionization directly to the 
continuum, is given by the FC factor between the 
ground level and the vibrational levels of the ionic state. 
The /3 values are expected to be independent of the 
vibrational quantum numbers of the ionic state. If 
photoelectron emission occurs through the resonance, 
on the other hand, one should take account of both the 
resonant and non-resonant components with the help of 
Fano's resonance theory [43]. Then the branching ratio 
is governed not only by the FC factor between the 
ground level and the vibrational levels of the ionic state 
but also by the FC factor between the vibrational level 
of the excited state and the vibrational levels of the ionic 
state. If the resonant enhancement is so strong that the 
non-resonant component is negligible (i.e., \q\ » 1 in 
Fano's resonance formula [43]), the branching ratio is 
dominated solely by the FC Factor between the vibra- 
tional level of the excited state and the vibrational levels 
of the ionic state. In this situation, the /3 -values are 
again independent of the vibrational quantum numbers 
of the ionic state. 
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Fig. 13. Plots of the branching ratio for the (320), (420), and (520) levels of the COJ ground state. 
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Fig. 15. Plots of the branching ratio for the (001), (101), and (201) levels of the COJ ground state. 



813 



Volume 106, Number 5, September-October 2001 

Journal of Research of the National Institute of Standards and Technology 



301 



0.05 




0.05 



0.04 



680 



700 



501 




720 



740 



760 



Wavelength [A] 



780 



800 



Fig. 16. Plots of the branching ratio for the (301), (401), and (501) levels of the COJ ground state. 
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Fig. 17. Plots of the branching ratio for the (211), (311), 411), and (511) levels of the COJ ground state. 
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In our previous publication [15], we made assump- 
tions mentioned above and the approximation of neglect- 
ing the contribution from the non-resonant ionization on 
the basis that the line shapes were in general Lorentzian. 
In terms of the Fano [43] theory, this means that the 
Igl-value is very large, indicating that the continuum 
component is very small. Then the FC factors for the 
transition between the vibrational level of the resonant 
state and the vibrational levels of the final ionic state 
were calculated. The approximation was made that the 
potential parameters for the resonant state are the same 
as those for the ionic state to which the resonance 
converges. Clearly this approximation should improve 
the higher up the Rydberg series the resonance is, 
and some qualitative agreement with this method of 
calculation was found. For the (0,0,0) branching ratio, 
the A«i and A« 2 (« = 4-7) TO series are calculated to 
have a slowly increasing branching ratio as one proceeds 
up the Rydberg series approaching the theoretical value. 
However, there was poor agreement for the branching 
ratios as one proceeds along the members of the main 
vibrational progression (n„,0,0), as noted in [15], and it 
was clear that considerable intensity moves from this 
progression into other vibrational levels. 



To take account of the intensity transfer from the 
main progression to the other vibrational modes as a 
first approximation, we have added the intensities of the 
nearest vibrational member above and below the (v ,0,0; 
v = 1-7) vibrational members for the same resonances 
as in the earlier paper. For example, referring to Table 1 , 
we added the values for (0,2,0) and (0,0,1) to the 
branching ratio for (1,0,0) and similarly for the other 
vibrational levels up to (5,0,0). Our labelling of the 
nearest levels changes for (6,0,0) and (7,0,0), but, as 
pointed out earlier, this labelling is not definitive given 
the resolution of our experiment. It was derived from the 
fitting process, but it would be impossible to confirm 
that intensity found in the (3,1,2) level could not in fact 
be attributed to the (5,2,0) level which was not included 
in our model. In the case of the (6,0,0) and (7,0,0) 
levels we added in the intensities from the vibrational 
levels just above and below as with the other vibrational 
levels. 

The results are shown in Fig. 21 as a bar chart, where 
each group of intensities corresponds to the resonance 
shown. The heights of the bars correspond to the 
branching ratios for (0,0,0) to (7,0,0), going from left to 
right in each group, with the additional contributions 
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added as described above. In Fig. 22, the theoretical 
values from the previous work are shown. In principle 
these should be the same for any given Rydberg series, 
and are labelled accordingly. Some similarities are 
immediately obvious: 

• For the first set of resonances, An a , n = 4 - 7, 
although theory puts (1,0,0) more intense than 
(0,0,0), the intensity dies away as one proceeds up 
the vibrational progression much as theory suggests. 

• For the next set, A«i, n = 4-7, a minimum at 
(2,0,0) and a maximum at (5,0,0) is more or less 
reproduced. For A7i, the (0,0,0) intensity is about 
right, and here the theory looks quite good. 

• For the next set, An 2 , n = 4 - 7, a low value of 
(1,0,0) is not reproduced, but the minimum at (3,0,0) 
is there, together with higher values of (5,0,0), 
(6,0,0) and (7,0,0) compared with the other reso- 

0.4 



0.35 



nances. A5 2 , also looks particularly good; (0,0,0) is 
a bit low, but the overall structure is clearly there. 

In the last set, An 3 , n = 4 - 7, two minima are pre- 
dicted at (2,0,0) and (5,0,0). This fits well with A5 3 , 
and A7 3 , with some suggestion of it at A4 3 . Again 
theory puts (1,0,0) as most intense, which is not the 
case. 

The general trend for intensity to move to higher 
vibrational members as one moves to higher 
vibrational excitations in the resonances is also 
reproduced. This is evident in Fig. 5; in the 
resonance regions vibrational structure extends, with 
considerable intensity, well into the regions of 
high vibrational energy, in contrast to the off- 
resonance spectra where almost all the intensity is 
concentrated in the first few members of the main 
vibrational progression. 



O 
05 

o 

LL 



0.3 



0.25 



0.2 



0.15 



0.1 



0.05 




Theory 



Theory 



Theory 



Theory 



Fig. 22. Histogram of the calculated Franck-Condon factors for the COJ ground state at positions in the absorption 
spectrum where Rydberg autoionizing resonance levels exist. The calculation could not distinguish between the various 
members of a particular series and hence the same distribution is attributed irrespective of principal quantum number. The 
resonances calculated are for the Tanaka-Ogawa series with vibrational levels of 0, 1 ,2,3 quanta of the symmetric stretch in 
the ground electronic state of CO!- 



820 



Volume 106, Number 5, September-October 2001 

Journal of Research of the National Institute of Standards and Technology 



It seems, therefore, that our FC model, which 
neglects direct ionization and takes account of the 
intensity transfer to the nearest different vibrational 
levels, can be applied in resonance regions for some 
cases. We do not at the present time have any general 
rules for quantifying the intensity transfer to the nearest 
vibrational levels for all the spectral regions. Some 
disagreements may be attributed to the neglect of direct 
ionization. 

An analysis on the basis of that carried out by Smith 
[20] for 2 , in which the Fano parameters are derived by 
fitting the experimental line shape to the Fano formula, 
includes the direct ionization component. Given our 
assumption above based on Lorentzian line shapes, and 
the fact that there are many overlapping lines, it was 
not thought worthwhile to follow this procedure here, 
although for the Hopfield series in the region between 
680 A and 720 A where the line shapes are not sym- 
metric it may be more appropriate. Even here, however, 
our experiment has insufficient resolution to separate 
overlapping structure clearly and values of the Fano 
parameters obtained by any fitting procedure are likely 
to be in some doubt. 

Turning our attention to the measured /3 -values, we 
note that the average value of /3, where we define 
"average" as an underlying value of /3 outside resonant 
structure over the whole wavelength range (i.e., the 
regions corresponding to direct ionization), is between 
and -0.5 in general for all the spectra shown. Even for 
the angular distribution measurements not shown here, 
because there was considerable scatter in the data, this 
still appeared to be the case. This would be consistent 
with, though not evidence for, validity of the Born- 
Oppenheimer approximation and the FC principle for 
direct ionization. The fact that the j3 value is negative 
indicates the presence of "parity unfavoured" transitions 
[41]; an analysis on the basis of the angular momentum 
theory proved partially successful for the oxygen 
molecule [42], where only one vibrational mode is 
present. Given the larger number of decay channels 
available for C0 2 , an analysis on this basis is unlikely to 
be conclusive. 

In the resonance regions, however, we notice break- 
down of the FC model applied as described earlier, 
since that model dictates that the /3 -values should 
not change for vibrational levels within the same 
progression. Note for example in Fig. 6 the markedly 
different response of /3, in the case of the two vibra- 
tional levels shown, to the B 3 (s) resonance at —751 A. 
The breakdown of the simple picture may be just due to 
neglect of the non-resonant component in the region of 
the resonances. The B 3 (s) peak is however the 
strongest among all the peaks in the total cross section 
spectrum, and it is hard to believe that the non-resonant 



component can contribute to this dramatic effect. The 
different response of the different vibrational levels 
within the same progression suggests breakdown of 
the Born-Oppenheimer approximation and/or the FC 
principle. 

The results presented in this work show that there is 
intensity in the odd bending modes, for example 
(vi,l,v 3 ), contrary to the selection rules. The way in 
which the branching ratios and /3 -values for the bending 
modes respond to the resonances varies dramatically: 
the strong response of the (0,1,0) /3 -value to the B3 (s) 
resonance in Fig. 9, and in Fig. 10 the response of the 
(2,1,0) branching ratio to the B5 (s) resonance, the 
(3,1,0) branching ratio to the A5 3 (TO) resonance, the 

(4.1.0) branching ratio to the A5i(TO) resonance and 
the (5,1,0) branching ratio to the A5 2 (TO) resonance. 
However, our experiment cannot resolve these levels 
from the corresponding (0,0,2), (1,0,2), (2,0,2) and 
(3,0,2) levels, so it is not clear the selection rules are in 
fact being broken. 

In general the response seems unpredictable; it does 
not appear to follow the generally accepted vibrational 
transition relative intensity expectations. For example, 
in Fig. 17 the branching ratios for the overtone modes 

(2.1.1) and (3,1,1) are enhanced at the B8 (s) resonance. 
The statistical errors in the corresponding values for /3 
were too large for these modes to show any structure; 
note that the branching ratio is down at the 4 % level or 
less. On the other hand, in Fig. 18 where the branching 
ratios are down at the 2 % level, it was possible to 
discern some structure in the /3 -parameter measure- 
ments; the ability to do this depended very much on the 
quality of the fit, discussed above. However no clear 
pattern emerges, and it is not at all obvious what the 
connection is between one particular branching ratio or 
j3 -value being enhanced, and the resonant structure. 
Given the uncertainties of the identifications of the 
weaker vibrational levels, the conclusion must be 
that weak modes such as («,1,0) are not necessarily 
enhanced by the resonances, a conclusion which is sup- 
ported by the results of a higher resolution experiment 
which will be the subject of a forthcoming publication. 

The detailed analysis of these data remains an out- 
standing task. The complicated nature of the inter- 
actions, and the fact that many of the resonances overlap 
or are incompletely resolved, makes theoretical 
modelling of these data very difficult. Nevertheless, one 
point emerges very clearly: it is necessary to make 
measurements of both the dynamical variables, the 
branching ratios and /3 parameters, in order to reveal 
weak structure not evident in the photoionization 
spectrum. The differing response of these two para- 
meters is particularly evident in Fig. 6 for the (0,0,0) and 
(1,0,0) members. Furthermore, for the weaker modes, it 



821 



Volume 106, Number 5, September-October 2001 

Journal of Research of the National Institute of Standards and Technology 



is necessary for the experiment to have sufficient 
sensitivity to measure branching ratios at the 1 % level, 
with superior resolution to that for the data shown here, 
in order to be definitive. 

The data described in this paper are available on the 
NIST web site at http://physics.nist.gov/co2arpes 

Acknowledgments 

The support of the staff at the Synchrotron Radiation 
Source at Daresbury Laboratory in the UK while 
making these measurements was much appreciated. 
This work was funded in part by the Science and 
Engineering Research Council in the UK. 

5. References 

[1] A. Dalgarno and J. L. Fox, Unimolecular and Bimolecular 

Ion-Molecule Reaction Dynamics, C. Y. Ng, T. Baer and I. 

Powis, eds., Wiley series in Ion Chemistry and Physics, Wiley, 

Chichester (1994). 
[2] C. A. Barth, C. W. Hord, J. B. Pearce, K. K. Kelly, G. P. 

Anderson, and A. I. Stewart, J. Geophys. Res. 76, 2213 (1971). 
[3] Y. Tanaka and M. Ogawa, Can. J. Phys. 40, 879 (1962). 
[4] Y. Tanaka, A. S. Jursa, and F. J. LeBlanc, J. Chem. Phys. 3, 

1199 (1960). 
[5] E. Lindholm, Ark. Fys. 40, 129 (1969). 
[6] C. Fridh, L. Asbrink, and E. Lindholm, Chem. Phys. 27, 169 

(1978). 
[7] H. J. Henning, Ann. Phys. (Leipzig) 1, 599 (1932). 
[8] K. E. McCulloh J. Chem. Phys. 59, 4250 (1973). 
[9] M.-J. Hubin-Franskin, J. Delwiche, P. Morin, M. Y. Adam, 

I. Nenner, and P. Roy, J. Chem. Phys. 8, 4246 (1984). 
[10] D. M. P. Holland and M. A. Hayes, Chem. Phys. 148, 241 

(1990). 
[11] T. Baer and P. M. Guyon, J. Chem Phys. 85, 4765 (1986). 
[12] R. J. Buenker, M. Honigmann, H.-P. Liebermann and M. 

Kimura, J. Chem. Phys. 113, 1046 (2000). This paper has exten- 
sive references on CO2 research using electron energy loss 

spectroscopy. 
[13] K. Takeshita, N. Shida, and E. Miyoshi, J. Chem Phys. 112, 

10838 (2000). 
[14] M. Johnson and J. Rostas, Mol. Phys. 85, 839 (1995). This paper 

has extensive references to the theoretical and experimental 

spectroscopic literature on C0 2 . 
[15] J. B. West, M. A. Hayes, J. B. West, M. R. F. Siggel, J. L. 

Dehmer, P. M. Dehmer, A. C. Parr and J. E. Hardis, J. Chem. 

Phys. 104, 3923 (1996). 
[16] A. C. Parr, P. M. Dehmer, J. L. Dehmer, K. Ueda, J. B .West, 

M. R. F. Siggel and M. A. Hayes, J. Chem. Phys. 100, 8768 

(1994). 
[17] A. C. Parr, D. L. Ederer, J. L. Dehmer and D. M. P Holland, J. 

Chem. Phys. 77, 111 (1982). 
[18] P. Roy, R. J. Bartlett, W. J. Trela, T. A. Ferrett, A. C. Parr, S. 

H. Southworth, J. E. Hardis, V. Schmidt and J. L. Dehmer, 

J. Chem Phys. 94, 949 (1991). 
[19] M.-J. Hubin-Franskin, J. Delwiche and P.-M. Guyon, Z Phys D 

Atoms, Molecules and Clusters 5, 203 (1987). 
[20] A. L. Smith, J. Quant. Spect. Rad. Trans. 10, 1129 (1970). 
[21] F. H. Mies, Phys. Rev. 175, 164 (1968). 
[22] J. H. D. Eland, J.Chem Phys. 72, 6015 (1980). 
[23] J. H. D. Eland, J. Chim. Phys. 77, 613 (1980). 



[24] D. M. P. Holland, J. B. West, A. A. MacDowell, I. H .Munro and 

A. G. Beckett, Nucl. Instrum. Methods B44, 233 (1989). 
[25] E. Harting and F. H. Read, Electrostatic Lenses (Elsevier, 

Amsterdam, 1976). 
[26] A. C. Parr, S. H. Southworth, J. L .Dehmer and D. M. P. Holland, 

NucI.Instr. Meth. 222, 221 (1984). 
[27] V. G. Horton, E. T. Arakawa, R. N. Hamm and M. W Williams, 

Appl. Opt. 8, 667 (1969). 
[28] G. V. Marr and J. B. West, At. Data and Nucl. Data Tables 18, 

497 (1976). 
[29] D. M. P. Holland, A. C. Parr, D. L. Ederer, J. L. Dehmer and 

J. B. West, Nucl. Instrum. Methods 195, 331 (1982). 
[30] J. B. West, Vacuum Ultraviolet Photoionization and Photodisso- 

ciation of Molecules and Clusters, C.-Y Ng, ed., World 

Scientific, Singapore (1991) Chapt. 8. 
[31] S. R. Stockbauer, B. E. Cole, D. L. Ederer, J. B. West, A. C. Parr 

and J. L. Dehmer, Phys. Rev. Lett. 43, 757 (1979). 
[32] P. M. Dehmer, J.Chem.Phys. 83, 24 (1985). 
[33] J. Berkowitz, Photoabsorption, Photoionization, and Photo- 
electron Spectroscopy (Academic Press, New York, 1979). 
[34] P. R. Bevington, Data Reduction and Error Analysis for the 

Physical Sciences (McGraw-Hill, New York, 1969). 
[35] D. Gauyacq, C. Larcher and J. Rostas, Can. J. Phys. 57, 1634 

(1979). 
[36] H. Veenhuizen, B. Wannberg, L. Mattsson, K. E. Norell, 

C. Nohre, L. Karlsson, and K. Siegbahn, J. Electr. Spectrosc. 

Rel. Phen. 41, 205 (1986). 
[37] L.-S. Wang, J. Reutt, Y T. Lee and D. A. Shirley, J. Electr. Spect. 

Rel. Phen. 47, 167 (1988). 
[38] G. Chambaud, W Gabriel, P. Rosmus and J. Rostas, J. Phys. 

Chem. 1992, 3285 (1992). 
[39] B. Taylor and C. Kuyatt, Guidelines for Evaluating and Express- 
ing the Uncertainty of NIST Measurement Results, NIST 

Technical Note 1297 (1993). 
[40] R. B. Cairns and J. A. R. Samson, J. Opt.Soc. Am. 56, 1568 

(1966). 
[41] D. Dill, Phys. Rev. A7, 1976 (1973). 
[42] A. C. Parr, J. B. West, G. V. Marr and J. L. Dehmer, J. Phys. B: 

At. Mol. Opt. Phys. 31, 5161 (1998). 
[43] U. Fano, Phys. Rev. 124, 1866 (1961). 

About the authors: Albert Parr is a physicist and is 
chief of the Optical Technology Division in the NIST 
Physics Laboratory. John West is a physicist and is head 
of the vacuum ultraviolet program at the Daresbury 
Laboratory synchrotron radiation facility. Michelle King 
is a chemist and was a postdoctoral assistant on this 
project at the Daresbury synchrotron. Kiyoshi Veda is a 
physicist from Tohoku University and participates 
in atomic and molecular programs at the Daresbury 
synchrotron. Patricia Dehmer was a chemist at Argonne 
National Laboratory when this work was performed and 
is now the director of the Office of Basic Energy 
Sciences for the Department of Energy. Joseph Dehmer 
was a physicist at Argonne National Laboratory when 
this work was performed and is now the director of the 
Physics Division at the National Science Foundation. 
The National Institute of Standards and Technology is 
an agency of the Technology Administration, U.S. 
Department of Commerce. 



822 



